Time ordering and counting statistics 



S. Bachmann G.M. Graf G.B. Lesovik ^ 

^ Theoretische Physik, ETH-Hdnggerberg, 8093 Zurich, Switzerland 
( b > L.D. Landau Institute for Theoretical Physics RAS, 117940 Moscow, Russia 

September 21, 2009 

Abstract 

The basic quantum mechanical relation between fluctuations of transported charge 
and current correlators is discussed. It is found that, as a rule, the correlators are to be 
time-ordered in an unusual way. Instances where the difference with the conventional 
ordering matters are illustrated by means of a simple scattering model. We apply the 
results to resolve a discrepancy concerning the third cumulant of charge transport across 
a tunnel junction. 

1 Introduction 

Transport in mesoscopic systems has been discussed using different approaches, based on cur- 
rent correlators [HE], [3], either time- or in/out-ordered, or on the statistics of the transferred 
charge [I] (and, related to the latter, on the precession of a spin coupled to current [5]). In 
this work we intend to address two main points. First, we describe the relation between the 
counting statistics and the time ordering of correlators. It turns out that the correct time 
ordering differs as a rule from the conventional one, T , and is given by the Matthews' T* - 
ordering [6]. Second, we will present a model of energy independent scattering, where that 
difference matters. Though the model is implicit in previous works [21 El CJ , its formalization 
allows to establish the equivalence between the in/out- and the T* -ordering of currents and 
hence between the two approaches mentioned at the beginning. 

An application of these findings is the clarification of a discrepancy between [S] and [7J 
concerning the third cumulant of charge transfer through a tunnel junction. In [7] the dis- 
crepancy was claimed to be entirely due to the difference existing between unordered and 
T -ordered correlators [SI HI- Our explanation is different, as we spell out momentarily. 

An often chosen framework rests on states \t, a) forming a basis and labelled by their time 
t of passage across a fiducial point and w.r.t. a reference dynamics, as well as by further 
quantum numbers a. To be precise, in such a scheme t is an observable which entails as its 
canonically conjugate operator a Hamiltonian whose spectrum similarly covers the real line. 



The description is therefore an effective one valid near the Fermi energy; indeed, from that 
perspective the energy spectrum appears unbounded. This remark being made, the scattering 
amplitude from to 1*2,0:2) is denoted by S a2ai (t 2 — ti). Equivalence between in/out- 

and T -ordering has been shown [TTJJ [7j under the assumption that the scattering matrix 
satisfies 

S(t) = 0, (t<0). (1.1) 

The assumption, there called causality, should be called strict causality, as it e.g. rules out the 
limiting, but simple case of instantaneous, S(t) oc 5(t), or equivalently, energy independent 
scattering. Actually, in the latter case, disagreement between the two orderings was found [HI 
[TT] . Our observation that T* rather than T matters explains this difference, in that only for 
strictly causal scattering the two orderings agree. Generally the difference takes the guise of 
a Schwinger term related to the infinite depth of the Fermi sea, at least for the lowest order 
cumulant where the difference matters, i.e. the third. That cumulant has been the object of 
experimental work [T2l [13], which may be read as a confirmation of the T* ordering in the 
above mentioned model. 

The plan of the paper is as follows. In Section [2] we discuss the general relation between the 
generating function of the moments of the tranferred charge and the time ordering of current 
correlators. It is unrelated to scattering and not necessarily placed within a scheme of second 
quantization. In Section [3] we introduce a model of energy independent scattering at the level 
of first quantization. We also discuss the equivalence between in/out- and T* -ordering. In 
Section H] we will promote the model to second quantization. The quantities of Section [21 as 
well as others, can then be computed explicitly and the relation found there illustrated on 
this example. 

2 The T* -ordering of current correlators 
2.1 The result and its context 

We consider a current carrying device as symbolically illustrated in Fig. 1. The generating 
function of counting statistics [U] is 

x(A,t) = (e^) = ^^L((Ag) fe ), (2.1) 

fe=0 

where ((AQ) k ) is the k-th moment of the charge transported during time t. Similarly, logx 
generates the cumulants (((AQ) fc )). Several quantum mechanical expression for x have been 
proposed, based in part on different measurement protocols. The first proposal [9] is 

X (\,t) = (exp(iXJI(t')dt')) (2.2) 
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Figure 1: A device connecting two leads. 

where I(t) is the current operator at the junction and at time t, and (•) has become the 
expectation in the initial quantum state of the system. The second one [I] reads, as recast by 

x (\,t) = ( e iA Q(*) e - iA «) = ( e iHt e iXQ e- im e- iXQ ) , (2.3) 

where Q(t) = exp(iHt)Q exp(— iHt) is the charge operator on the right of the junction and 
at time t, and H is the Hamiltonian. The definition is appropriate to the situation where 
the initial state is an eigenstate of Q and the charge is measured after time t. In fact, if the 
initial eigenvalue is q, then x(^M) = (e lXQ ^}e~ lXq just describes the statistics of q' — q, where 
q' is the (random) outcome of the measurement of Q(t). The quantity q' — q is identified with 
the transported charge AQ. We mention in passing that there is a generalization [15] of this 
definition to the case where the assumption on the initial state does not apply and q is the 
outcome of the measurement of Q(0) = Q. 

A further proposal [TBI E] is based on observing a spin coupled to current 

X (A, t) = (e^(-V2)* e -iH(A/2)tj = ( e -i±Q e iHt e iXQ e -iHt e -i^ ? ^ 

where H(\) = e lX ®He~ lX ® . This approach does not require the initial state to be an eigenstate 
of Q, but if it does, then it agrees with (I2.3p . A related proposal was put forward in [S] (with 
similar expressions found in [171 ITU] ). It is given by the Keldysh time-ordered expression 

X (X,t) = (f exp (i(A/2) / I(t')dt')T exp (i(A/2) / I{t')dt')) , (2.5) 

o o 

where T denotes the usual time-ordering and T the one in the opposite direction. For later 
purposes we recall that the time ordering is supposed to occur inside the integrals once the 
exponential is expanded in powers of A. 

Eqs. (12 .4j) and ( 12. 5ft may differ in some applications, as it is the case in the very simple 
situation of independent fermions passing with fixed transparency T across a scatterer biased 
by V . In fact it was found in [5] and [S] that the third cumulants are, respectively, 

(((AQ) 3 )) = T(l - T)(l - 2T)^ , (((AQ) 3 }} = -2T 2 (1 - T)^ . (2.6) 

Z7T Z7T 

In [7J the first answer has been found for the second approach, i.e. for (12. 5ft . by trading the 
Keldysh ordering for in/out-ordering. This step in turn relies on the strict causality of the 
scattering matrix, which however is not consistent with an energy-independent transparency. 
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The results of this article explain the discrepancies. We will consider (12. 3p for simplicity, 
though similar considerations could be made for (12.41) . As mentioned before, the two defini- 
tions agree in some cases, including the model underlying (12.61) . We want to point out that 
some care is to be exercised when expressing the moments in terms of time-ordered current 
correlation functions. The result is 

<(AQ) fc >= f t d k t{T*(l(t 1 )---I(t k ))), (2.7) 
Jo 

where d k t = dti---dt k , I(t) = dQ(t)/dt, Q(t) = exp(iHt)Q exp(— iHt), and T* means that 
the derivative has to be taken after the time ordering: 

T* (I( tl ) ■ ■ ■ I(t n )) := A • • • A T (Q( tl ) ■ ■ ■ Q(t n )) . (2.8) 

Hence, eq. Q may be summarized as x(A,t) = (T* exp (i(A/2) I{t')dt'))- by (TJS]), this 
amounts to keeping the time ordering outside of the integrals of the exponential series. Sim- 
ilarly, (12 .4p is equivalent to (12.51) once stars are added to the time orderings. 

It is known P, [18] and will be shown below that, as a rule, (12. 8ft differs from T (/(ti) • • • 
I(t n )) by contact terms supported at coinciding times but agrees with it if [Q, I] = 0. 



2.2 The derivation 

In this section we show (12.71) and obtain an expansion of (12.81) in contact terms. Any arising 
operators will be assumed to be properly defined and that will be implicitely verified in later 
sections in the context of applications. 
Either side of 

( (e^V^e-iH*) e -iAQ ) = ( e ltft ( e iAQ e -iift e -iAQ) ) (2 g) 

is an expression for x(^,t) an d may serve as a starting point for the derivation of (12. 7\\ . We 
will follow both avenues, because the l.h.s. leads to a simple proof, while the r.h.s. has often 
been considered in the literature [5"| \T§\. 
We begin with the l.h.s., which yields 

x (A,t) = (e iA «We" iAQ ) = (Te iA(Q(t) - Q) ) (2.10) 

and hence 

<(AQ) fc > = (T (Q(t) - Q) k ) = (T ((Q(t0 - Q) ■ ■ ■ {Q{t k ) ~ Q)))\ tl= ... =tk=t ■ (2.H) 

We note that the expression in brackets vanishes if t j = for some i — 1, . . . , n. By repeated 
use of the fundamental theorem of calculus we obtain 

<(AQ) fc > = f dtx A (T ((Q( tl ) - Q) ■ ■ ■ (Q(t k ) - Q)))\ t2= ^ =tk=t 
= f dH A . . . |_ (T ((Q( tl ) -Q)-.. (Q(t k ) - Q))> . 
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Expanding the correlator in Q(tj) and —Q, the resulting second term is independent of tj and 
does not contribute to the derivative. This proves (12.71) . 
Let us now turn to the r.h.s. of (12. 9ft : 

X (X,t) = (e iHt e- iH ^) , (2.12) 

where 

H(X) = e iXQ He~ iXQ = jr ^ad^(H) (2.13) 

j=0 J ' 

with multiple commutators defined by a,d° A (B) = B, ad^(_B) = [A, ad^ _1 (.B)], (j > 1). There- 
fore, e 1 ^* e - 1 ^( A ) f i s the propagator in the interaction picture for H with interaction 

W = H(X) - H = ijr ^ad^I) • (2-14) 

It is expressed by the Dyson expansion 

00 (—\) n ft 

e im e -iH(x)t = L_L / d n tT(W{t 1 )---W{t n )) (2.15) 

with W(t) = exp(iHt)W exp(-iHt) . The term of order X k in the expansion is obtained by 
picking the power A Ji from W(ti) through (12.141) . in such a way that YH=\3i = ^ : 

((AQ) k }= Yl ^ ■ , r^t(T(ad^- 1) (/)(t 1 )---adg"- 1) (J)(t n ))). (2.16) 

The domain of integration is invariant under permutations of the times and the integrand 
under those of its factors. Thus, n-tuples . . . ,j n ) differing only by order may be binned 
together. To this end, let 

nj = m\ji=j} (2-17) 



be the number of times the power A J has been picked. Each bin then consists of 



nil ■■■n k \ 

tuples, since permutations among equal j;'s do not generate new tuples. Thus, 

-t 



(2.18) 



(fH.na,...) (JL=l^!J(JI i= lj! Uj ) J ° 



(2.19) 

where n = Ylj n j an< ^ the iz' s satisfy (12.171) . The products are finite because of the condition 

The derivation just given establishes the equality between ( 12.191) and (12.71) . In the following 
we give an independent one by expanding ( 12 .7\\ in contact terms. 
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Proposition 1. 



T* (1(h) ■ ■ ■ I(t k )) = A • • ■ A T (Qfo) ■ • • QftO) 



E T (lI ad Q CM ( J )(^)^)' ( 2 - 2 °) 



where the sum runs over all partitions P of {1, ... ,k} into nonempty disjoint subsets C. 
By 5c we understand a product of 5 -functions collapsing the times ti, (i G C) to a single 
time tc = ti. More precisely, for C consisting of n(C) elements, C = {i\, . . . , z n (c*)}; we se t 
$c = Tiffi' 1 5(tij - tij+i); in particular 5 C = 1 for n(C) = 1. 

Remark. The partition into single-element clusters contributes T (I(ti) ■ ■ • I(tk)). 
Proof. We claim the slightly more general statement for < I < k: 

A • • • J-T (Qih) ■ .-Q(t k )) = T (( II ^ C) -\mc)Sc) ■ Q{ti + x) ■ ■•<?(**)) , (2.21) 

1 1 PGP, ceP 

where the sum now runs over V\ instead of Vk- This is trivially true for I = and identical 
to (I2.20p for I = k. To run the induction in I, we note that for operators A{ associated with 
times ti, 

r) T 

—T{A 1 ---A r ) = T(A 1 ---A j ---A r )+ S(t j -t i )T([A j ,A i ]A 1 ---A i ---A r --A r ), 

^ i=l, i^j 

(2.22) 

where " denotes omission. Using this with j = I + 1 on ( 12. 2 1 1) we observe that 5(ti + i — 
ti)[Q(ti+i), Q(ti)] — 0. The other commutators so generated are 

S(t l+1 -t c )[Q(t l+1 ), a d n Q (c) -\l)(t c )] = 5{t l+l -t c )^ c \l){t c ) , (2.23) 



so that 

d d 



T(Q(t 1 )'--Q(t k )) 



dt l+1 <9ti 

= E T ((II ^f ] ~ l mc)5c) ■ I(t l+l )Q(t l+2 ) ■ ..Q(t k ] 



PeVi cgp 



EE T Ud n Q {c '\l)(t c ,)5 c ,5(t l+1 -t c ) ■ ( J] a d n Q W-\l)(tc)6c) 



PeVi C'eP CeP 

■Q{ti+2)---Q(K 

This agrees with (12.211) for /+1 in place of I. In fact, partitions of {1, . . . , 1+1} are distinguished 
by whether the cluster C containing / + lisC = {/+l} (first line) or C = C U {/ + 1} with 
C e P eVi (second line); in this case, n(C) - 1 = n(C). □ 
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We can now verify the stated equality. After inserting (12.201) in (12.71) . partitions P G Vk 
can be binned according to the numbers rij, (j = 1,2, . . .) of their clusters of size j (hence 
Y2jJ n j = Each bin consists of 

fc! 

mMWJ^ (2 ' 24) 

partitions. In fact, permutations of k elements permuting or preserving clusters do not gener- 
ate new partitions. Since in (12.201) . J2 c n(C) = k and ^ % = n, the equality between (12.191) 
and (12 .TP is established. 



3 The instantaneous scattering model 
3.1 The model 

Transport of charge across a scatterer at zero temperature and at a small bias ought to be 
determined fully by the Fermi velocity (equal to 1 in suitable units) and the scattering matrix 
at the Fermi energy. The Fermi sea is, effectively, infinitely deep. The model we are about to 
define is minimal in the sense that it is fully determined by that matrix. The left and right 
reservoirs are represented by two infinite one-dimensional leads, which are however chiral. 
The scattering between leads occurs instantaneously at a single point. The model may be 
seen as describing a quantum point contact. 

Let H = L 2 (R) © L 2 (R) = L 2 (R;C 2 ) be the single-particle Hilbert space of the model. 
The charge operator Q is the projection onto the right lead, 

«-(!;)• ™ 

the second quantization of which will later represent the number of particles there. In ab- 
sence of interaction, states i/j G Ti evolve in time according to i])(x — t). The corresponding 
Hamiltonian H Q is linear in the momentum: 

where p = —id/dx and the derivative is taken in the sense of distributions. Scattering between 
leads is determined by the matrix 

S=(\l). (3.3) 



We will specify the resulting Hamiltonian [20] at the end of this section. First however, we 
define the unitary group generated by it, because it is simpler and it is all which is needed in 
the rest of the paper. It is 

(TT(A,IAM _j ^~t) + (S- 1M0 < x < t)^{x -t), (t > 0) , 

1 l JWl J ~ 1 i>(x-t) + (S*-l)9(t<x<0W(x-t), (t<0), [ ' 
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and is motivated by the idea that the part of the freely evolved wave function ip(x — t) which, 
say for times t > 0, has crossed the scatterer at x = between times and t gets replaced by 
Sip(x — t). One readily verifies that U(t) is a strongly continuous 1-parameter group. 
Letting U (t) act on the charge operator, we get 

nm - TT(A*DTT(A - S Q ^~ X ~ *) + 9 W + S *Q Sd ( x + *)*(-*) »(*>«), 

y W _ (7 W ^(7 [t) - | + e{x + f)) + S q S * 9 (_ x _ t)9{x) ) (t < o) , W 

or equivalently 

Q(t) = (Q6(-x -t) + S*QS6(x + t))6(-x) + (SQS*6(-x -t) + Q6(x + t))9{x) . (3.6) 
It yields the current 

lit) = = ((S*QS - Q)9(-x) + (Q - SQS*)8(x))5(x + t) , (3.7) 

which is well-defined distribution in t. 

We now come to the description of the generator H, a self-adjoint operator by Stone's 
theorem. To this end, let 

V + = G L 2 (R + ; C 2 ) | tfj' + G L 2 (R + ; C 2 )} , (3.8) 

where the derivative is that of a distribution on M + = (0, oo), i.e. away from the origin. 
Any function ip G V + is continuous up to the boundary x — 0, and let ^+(0) be its boundary 
value. Any function in L 2 (M + ; C 2 ) may be seen as a distribution on all of K. In this sense, 
the derivative of ip+ £ *D+ is 

^ + = 1 J/ + (x)+^ + (0)S(x). (3.9) 
We may similarly define G X>_ based on K_ = (— oo, 0). Then, 

-^iP_ = ijj'Jx) - ip(0)6(x) . (3.10) 
ax 

Given ip G L 2 (R; C 2 ), let ip± be their restrictions to the half-lines R±. 
Proposition 2. The generator H ofU(t) has domain 

D(H) = {ip £ H \ ip± £ V±, ^+(0) = SV>-(0)} (3.11) 

and is given by 

(Hrl>)(x) = -i(il>' + (x)+il>'_(x)). (3.12) 

We recall that D{H) consists of those states ip G Ti. for which the limit lim^o t^iJJit) — 
1)^ = —iHip exists, thereby defining Hip. 
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Proof. Let first t > 0. Then 

(U(t)ijj)(x) = i/j+(x -t)+ ifj^x -t) + (S- l)9(x)ip_(x - t) . (3.13) 
Using (13~91I3TT01) . we find 

= -W4 +<!>'-) + (SV'-(O) - i>+(0))6 (3.15) 

in the sense of distributions if ip± G T>±. Actually, this last condition is implied by the 
requirement tp G D(H). In fact, by using t _1 (0, (U(t) — l)ip) = t _1 ((£/(t) — 1)0, ■0) on a test 
function <fi supported away from x = 0, we obtain in the limit — i(0, i7"0) = + <p'_, ip) and 
hence ^ G L 2 (M ± , C 2 ). At that point ip G further implies Sip-(0) = ip+(0). Moreover, 

the convergence (13.151) is attained in 7i. The same conclusion in reached for t < 0, whence 
(JSISD is -iHvj). □ 
More casually, the Hamiltonian can also be written as 

H = H + i{S - 1) ■ <L = H - i{S* -l)-6 + = H + 2i^^ • 5 , (3.16) 

with quadratic forms 5± on V + n V- defined as (5±ip)(x) := 5(x)^(0±) and 5 = (5+ + 5J)/2. 
The first expression in (13.161) corresponds to (13.141) and the remaining ones follow from the 
boundary condition. Eq. (I3.16P is a variant of one found in [2Tj . 

We note in passing that ( 13. 7ft may alternatively be understood as a quadratic form on 
D(H) for fixed £. The expression is unambiguous even at t — 0, because (S*QS — Q)5- and 
(Q — S'QS'*)^ agree as a result of the boundary condition. 



3.2 Comparison of time orderings 

The T* -ordering of currents is identical to the in/out-ordering, T, whose definition [TU] 
we shall recall momentarily. For the time being, we establish this fact at the level of first 
quantization, but we shall show in Section 4 that it persists under second quantization. The 
current ( 13.71) may be split, I(t) = I+(t) + I-(t), into outgoing and incoming parts, 

/+(*) = (S*QS6(-x) + Q0(x))6(x + t), 

I_(t) = -(Q6(-x) + SQS*6(x))5(x + t). 1 ' ] 

In particular, [I±(t) , I±(s)] = 0. 
Proposition 3. We have 

T* (1(h) ■ ■ ■ I(t k )) = f (1(h) ■ ■ ■ 1(h)) , (3.18) 

where the ordering T places any I-(t) to the right of any I+(t'), regardless oft^t' (the order 
of currents of the same type is irrelevant), and extends by linearity to I(t) = I+(t) + I-(t). 
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Before giving the proof, let us make an observation. The proof of (12.201) was based on 

(I(f)Q(*)) = T (I(t)I(s)) + 5(t - s)[Q(s), 1(f)] , (3.19) 
which is a sensible equation in the context of our model. The formal analogue 

^T (I+(t)Q(s)) = T (I + (t)I(s)) + 5(t - s)[Q(s), I + (t)\ (3.20) 

is meaningless on its r.h.s., and should not be used. These claims are based on the commu- 
tators 

[Q(s),I+(t)] = e(t-s)5(x + t)([Q,S*QS]6(-x) + [SQS*,Q]6(x)), 

[Q(s), !_(*)] = e(s-t)8(x + t)([Q,S*QS}9(-x) + [SQS*,Q}9(x)), (3.21) 

[I(s),I+(t)\ = ~6(t-s)5(x + t)([Q 1 S*QS}6(-x) + [SQS*,Q}e(x)) , 

which can be obtained from (13.61 13 . 1 71) . The first commutator is discontinuous at t — s = 0, 
hence its multiplication with 5(t — s) in (I3.20p is ambiguous. Also the first term there is 
seen to exhibit an ambiguity proportional to 5(t — s). By contrast, the sum of the first two 
commutators, 

[Q(s), I(t)] = 5(x + t)([Q, S*QS]6(-x) + [SQS\ Q]6(x)) , (3.22) 
is independent of s, whence the r.h.s. of (13.191) is well-defined. 

Proof. The inductively stable generalization of (I3.18P is 

— ■ ■ ■ — T (Qitt) ■ ■ ■ Q(t k )) = f (1(1,) ■ ■ ■ Ifa); Q(t l+1 ) ■ ■ ■ Q(t k )) , (0 < I < k) , (3.23) 

where, on the r.h.s., the Q's are time-ordered and the i_ (resp. i+) are placed to their 
right (resp. left). Upon differentiating the r.h.s. w.r.t. t t+1 , eq. (I2.22p generates equal time 
commutators only among Q's, which vanish. Thus, 

^-f (/(tO • • • J(t,); Q(t l+1 )Q(t l+2 ) ■ ■ ■ Q(t k )) = f (J(t0 • ■ ■ /(*,); I(t l+1 )Q(t l+2 ) ■ ■ ■ Q(t k )) , 

(3.24) 

indicating that I(t/+i) is still subject to time ordering relatively to Q(t[ +2 ) ■ ■ -Q(t k ). After 
splitting it, I(t [+ i) = I + (t [+ i) + J_(t i+1 ), the incoming part I^(ti +1 ) commutes with the earlier 
Q's on its right by (I3.2ip ; similarly, J + (t i+1 ) with those on its left. The result is (I3.23P with 
/ + 1 instead of /. □ 



4 Second quantization 
4.1 Binomial statistics 

For fermionic many-body systems consisting of independent particles, the generating func- 
tion (12.31) has been computed by Levitov and Lesovik [1] as 

X (X, t) = ((1 - T) + e lX T) Vt/2n • (1 + o(t)) , (t - oo) , (4.1) 

10 



where T = |t| 2 is the transmission probability and V is the bias across the scatterer. The 
result describes a binomial distribution with N = tV/27r attempts. In particular, it yields the 
third cumulant 

<<AQ% = ^T(1-T)(1-2T). (4.2) 

Before presenting a derivation of (14.11) among others [HE2], let us make a digression on second 
quantization. In vague terms, the second quantization A of a single-particle operator A is the 
sum of its contributions from all particles i, A = — {^i)p)i though with expectation 

value subtracted. The notion can be formalized for infinitely many particles in the frame of 
the GNS space of the multi-particle state p, i.e. a Hilbert space containing p and its local 
perturbations as vectors. We shall do for p being a quasi-free fermionic state, determined by 
a single-particle density matrix of the same name, 1 > p = p* > 0. In the special case of 
a pure many-particle state, i.e. for p = p 2 , A admits a second quantization A acting on the 
GNS space iff 

[A,p]El 2 , (4.3) 

where X 2 denotes the Hilbert-Schmidt operators. The exponentiated version of this criterion 
is: A unitary operator U admits an implementation U on the GNS space iff [U, p] G I 2 - The 
notation is slightly abusive as the construction differs from A: The propagator U = exp(-iHt) 
is promoted to U = exp(— iHt). With these preliminaries, the generating function (12.31) reads 

X (X) = (u*e iX ®Ue-' lX ^ , (4.4) 

where, in line with the hypothesis made there, we assume 

[p,Q] = 0. (4.5) 

The generating function (14.41) can be expressed in terms of the first quantized operators 
through the Levitov-Lesovik (infinite) determinant [I]. In [22] a particle-hole symmetric 
variant of the determinant was given a mathematical foundation on the basis of the GNS 
space. We recall the result, in the case of a pure state. Set Au := U*AU. 

Proposition 4. Let 

[U,p]eZi, (4.6) 
where T\ (c X 2 ) denotes the trace class operators. Then 

X(X) = det (1 + (e" iA - l)Q uPuP ' + (e iA - l)Q uP ' uP ) , (4.7) 

where p' = 1 — p is the density matrix of holes. The determinant is Fredholm. 

In the context of (14.11) and of the instantaneous scattering model (13.161) . the initial density 
matrix is the projection 
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with pi = 6(ni — p), (i — L, R) representing two Fermi seas biased by V = pi — p R . While 
assumption (14.51) holds true, (14.61) unfortunately does not. The failure is however of minor 
importance. It can be traced back to the discontinuities of (U(t)ip)(x) at x = and x = t 
introduced by U, as seen in (13 .4p . Technically, the hypothesis could be met by making 
the scattering matrix time-dependent with S(r) = S for most of the relevant time interval 
r G (0, t), but = 1 near its ends (see [23J for a similar model). The particular choice of rounding 
would affect the determinant only to sub-leading order in t (cf. a logt contribution to x(A) 
in [5]). We neglect such contributions from the result, which in view of 9Ay9 = 9S*AS9 for 
9 = 9(— t < x < 0) becomes 

X (A, t) = det n{t) (1 + (e~ iA - l)S*QpSp> + (e iA - l)S*Qp'Sp) , (4.9) 

where H(t) = L 2 ([-t, 0]; C 2 ). The determinant is now of a form considered by Kac [21] and 
Akhiezer [25], extended to the matrix case [26]: Consider a translation invariant operator A 
on L 2 (M; C n ) with Fourier multiplier A(k). For its truncation to [—t, 0], one has 

t f°° 

logdet w(t) (l + A) — ► — / dk logdet(l + A(k)) + o(t) , (t -»• oo) , (4.10) 

where the determinant on the r.h.s is that of an n x n matrix. We apply the result to (14.91) 
with n = 2 and 

STQpSpt = p RP ' L ( _ T t I ) , (4.11) 

as well as with p and p' interchanged. For p L > p R we have Prp' l = and p' r Pl = G{^r < 
p < pi). We so obtain 

log X (A, t) = ^ log(l + (e iA - 1)T) + o(t) , (t - oo) , (4.12) 

as announced. 



4.2 The third cumulant 

In this section we apply the expansion in contact terms, eq. ( 12.201) . in order to compute the 
third cumulant of the transported charge. The currents will have to be understood as second 
quantized operators, which we emphasize by adding a hat. Eq. (12. 7ft reads 

«aq 3 » p = f^t((X(T(h)T(t2)T{t a )))) p (4.13) 

J o 

and the result to be re-derived is (14. 2p . It agrees with the binomial statistics as derived again 
in the previous section and the result [7] based on the T time ordering; on the other hand, it 
should be contrasted with 

^ dH((T{l( tl )I(t 2 )I(t 3 )))) p = --2T\l-T) (4.14) 
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which is the result obtained for ((AQ 3 )) p in other schemes [9J [8] - The difference between ( 14.131) 
and (I4.14p is accounted for by contact terms. Indeed, in eq. (I2.20p with k = 3, the partitions 
P are (1)(2)(3), (12)(3), (23)(1), (31)(2), (123), which yields 

((AQ 3 » p = [ t dH((TT l T 2 T 3 )) p + 3 f dH((Th[Q 2 ,T 2 ])) p + f dt x ll[Qi,[QiM)) P , (^) 
Jo Jo Jo 

where Qi = Q(U), ij = I(tA. Besides of (I4.14p . we will show that the two further integrals 
are 0, resp. {Vt/2ii) ■ T(l — T), resulting in 

((AQ 3 » P = ^(-2T 2 (1 — T) + + T(l — T)) = ^T(l - T)(l - 2T) . (4.16) 

It should also be remarked in passing that [Q, p] = 0, whence the hypothesis about the initial 
state, discussed after eq. ( 12. 3D and underlying eq. ( 14.131) . is satisfied. 

Before going into the proper proof of these claims, let us recall the following elementary 
rules [27], which reduce computations of correlators and commutators of second quantized 
operators to the level of first quantization: 

(A) P = , 

{{AB)) P = (AB) p = tr(pAp'Bp) , (4.17) 
((ABC)) P = tr(pAp'Bp'Cp) - tr(pAp'CpBp) , 

with p' = 1 — p. The r.h.s. of (I4.17P may be cast as tv(pAp'Bp) = tr([p, A]p'[B, p\), which is 
finite for A, B enjoying (14.31) . Eq. (I4.17P implies 

[A, B] = pCB] + (ti(pAp'Bp) - tr(p'ApBp')) ■ 1 = \A\B] + S(A, B) ■ 1 , (4.18) 
where the last term is known as a Schwinger term. It implies 



A(t) = A(t)+ifdt'S(H,A(t'))l. (4.19) 
Jo 

In our case, the current I(t) = dQ(t)/dt is a distribution in t. In order to give rise to 
an operator on which the above may be applied, a test function is required, which will 
however remain implicit. This being said, the contact terms in (14. 1511 are, by ( 14.181 14.19p . 

{{Th[Q 2 X])) P = ((Th[ChM)) P and (([4, [QiMh = (([QiAQ^h]])) P - They are of impor- 
tance since the commutator 

[Q(t), I(t)\ = ([Q, S*QS}6(-x) + [SQS*, Q}6(x))5(x + t) , (4.20) 

resulting from ( 13.221) . does not vanish. 

The computation of the integrals ( I4.15P will make repeated use of the following expressions. 
The translation invariant density matrices p and p 1 have integral kernels 

p(x, y) = ■ D(x — y) , 

p'(x, y) = ■ D(x — y) 

HK ,yj 2vriy-x-i0 v yj 
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with 



Moreover, 



S*QS = ( ^ f_ r V [Q, S*QS] = ( ° fi Q * ) . (4.22) 



We compute the first integrand (14. 15ft by temporarily dropping the time ordering. 

(( W 3 » P = tr( W/ 3 p) - txihp'hphp) (4.23) 
1 tr(A(t 2 - t!)A(t 3 - t 2 )A{h - t 3 )) - tr(A(t 3 - t x )A{t 2 - t 3 )A(t 1 - t 2 )) 



(2vri)3 (t, -t 2 - iO) (t 2 - t 3 - iO) fa - t 3 - iO) 



(4.24) 



with 



A(U - tj) = (S*QS - Q)D(U - tj) = f ^ ^ ) Dfa - , (4.25) 
where we used (13.71) . That results in 

//TTTW - ~^ 2 (l-T) sin Vih - t 2 ) + sin F (t 2 - t 3 ) + sin ^ fa - tj) 

{{hhh))p ~ (tl _ * 2 _ i0 )(t 2 - 1 3 - io)(tx - * 3 - io> ' (426) 

where V = — [j,r. In fact, the first trace equals 

2iT 2 (l - T) sin V(t 3 - h) + (cyclic) (4.27) 

by the following argument. It may be expressed as ■ k=L R AijAjkA^. The word ijk = LLL 
involves a single diagonal entry of (14.251) with overall compensating phases. Its contribution, 
T 3 , cancels an opposite contribution from the similar term RRR. For LRR, we find 

t '| e iMi{(*2-ti) . f_yWMfl(*3-*2) . ^ e i/4t(*i-*3) _ _y 2 (q — y) e _n/ (*3-*i) ? (4.28) 

which can be combined with T 2 (l - T)e lV ^-^ from RLL. The remaining words provide 
the cyclic permutations. To conclude, it suffices to note that the second trace is obtained by 
exchanging 2 <-> 3, which induces a change of sign in all the exponents. 

We then note that (14.261) is locally integrable and we drop the regularizations. Eq. (14.141) 
then follows from 

dh Sin V{h ~ h) + Sin J {h ~ * 3) + ™ V{h - h) = 2**Vt + o(t) , (f-oo). (4.29) 
{ti-t 2 ){t 2 -t 3 ){t 1 -t 3 ) W ' V ' V 7 

Actually, the time ordering should have been reinstated, but since the expression is permu- 
tation symmetric that is superfluous. Eq. (14.291) can be derived by reinstating the regular- 
izations — iO; then the integral can be broken into three terms, the first two of which vanish 
in the limit of large t. In fact, the first term has poles at t 3 — t 2 — iO, t\ — iO which do 
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not pinch the real axis. The same applies to the second and to the ^-integration. Using 
/ dt 2 {t l -t 2 - iO)" 1 • (t 2 - t 3 - iO)- 1 = 2m(t 1 -t 3 - i0) _1 , one is left with 

2vri / dhdt 3 - \ 3 1J = t ■ 2m / dx + o(t) = 2n 2 Vt + o (t) ; 4.30 

in fact, the odd part of (x + i0)~ 2 is ((x + i0)~ 2 - (x - i0)~ 2 ) /2 = iti6'(x). 
We next consider the middle integral in (14. 15ft . We find 

«£[&, / 2 ]»p = tr (V[Q 2 , J 2 ]p) = * _ i _ ■ tr(A(t a - *i)[Q, FQSlDfa - t 2 )) 

sinV^i - 1 2 ) , 



(2vri) 2 (*! - t 2 - iO) 2 



and the same result for (([Q2, h]h)) P except for +i0 instead of — iO. Thus, ((T7i[Q 2 , 1 2 ])) p is 
odd in t\ <-> t 2 , and the integral vanishes. 

Finally the last integrand in 04.15p . being the expectation value of a commutator, reduces 
to its Schwinger term 

[Qt, h]])) p = ti( P Q 1 p'[Q 1 , h]p) - tr(p'Qip[Qi, h]p') ■ (4.31) 
Using (jODl 133]) it equals 

L dy L + tl + w - { S+ u-w ) ^ s ' QsD{t > + »>»■ s ' qs w-'> - »» 



(27Ti) ! 



\t(1-T) / dj/smy^ + tO / 



(27ri) 2 v J J_ tl vy iy V(2/ + ^i+iO) 2 (2/ + ti-iO) 2 

where we noted that additional contributions from y < —ti and y > vanish because of 
tr(QD(ti + y)[Q, S*QS]D(—ti — yj) = 0. After changing the integration variable to x — y+h, 
both factors are odd, whence the integral equals 

\ C dx 8in Vx (WTW ~ ¥^w) = - %w + 0(1) ' (432) 

as in (14.301) . We conclude that 

J dti [Qi,h]])) P = T(l - T)^- , (4.33) 

as claimed. 



4.3 A strictly causal scattering process 

We may trade the instantaneous scattering process in use in the previous section with a 
strictly causal one, all while remaining within the model of Section [3j That is achieved by 
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simply forfeiting a piece of length I > of the leads in favor of the scatterer. So reinterpreted, 
the scattering process lasts 21 and the contact terms should disappear. Nevertheless the 
same result for the third cumulant will be obtained, though only the first terms on the r.h.s. 
of ( 14.15P will contribute. 

In physical terms we place the detector a distance / > away from the scatterer; in 
mathematical terms we replace Q by its regularization Qi = Q9(\x\ > I) and hence / by 
li = i[H,Qi] = Q[5(x — I) — 5(x + I)}. It ought to be noted that [Qi,p\ ^ 0, so that the 
physical appropriateness of (12.31) . and hence of (14.131) . could be questioned. Nevertheless, we 
may just view I > as a regulator before taking the limit I — > 0. However, even this is 
troublesome, at least in this form, because [Qi,p] T 2 makes Qi undefined. That can be 
remedied by smoothing the step function 6(\x\ > I) on a length < /. With these preliminaries 
taken, the contact terms in (14.151) vanish as expected. This follows from [Qi,Ii\ oc 1 and in 
turn from [Qi,Ii] = 0, while the corresponding Schwinger term no longer vanishes. On the 
other hand, the current correlator ((T I11I12I13}) in (14.151) has a well-defined limit when the 
smoothing is removed. Omitting the time ordering, it is for times tj > I [11] 

2 sin V{ti - t 2 ) + sin V(t 2 - t 3 ) + sin V(t 3 - tx) 



((I ll I l2 Iis)) P = T(l-T)(l-2T)- 
— T(l — T) 



(2vr)3 (ti-t 2 )(t2-t3)(*i-*3) 
2 / s\nV{t 1 -t 2 ) 



(2tt) 3 V (h ~ h)(h -t 3 -2l- iO) (h -h-2l- iO) 

siny(t 2 - h) 

(h -t 2 + 2l- i0)(t 2 - * 3 )(<i -h + 2l- iO) 

sin V(t 3 — ti) 
(ii - t 2 - 21 - i0)(t 2 -t 3 + 2l- i0)(ti - h) 



(4.34) 



= 1 + 11. 

We remark that the expression reduces to (14.261) in the limit / — > 0; the integral of the latter is 
(14.141) and should be distinguished from the limit of the integral. We add that time-unordered 



correlators, suitably symmetrized, can be measured by means of detectors discussed in [28J . 
The derivation of (14.341) begins as in (I4.23p . but with currents 

I(t) = S*QS5(x -l + t)- Q5(x + l + t) = /+(*) + /_(*) (4.35) 

split into outgoing and incoming parts. The terms involving only outgoing parts yield I. 
Indeed, their contribution has the structure (14.241 14.251) . but with S*QS — Q replaced by 
S*QS. In the description used there in reference to the first trace (14.241) . the words LLL and 
RRR contribute T 3 and (1 — T) 3 , respectively; LRR and RLL contribute T(l — T) 2 e~ iV '^ 3 "* l) 
and T 2 (l - T)e iV/( * 3 "* 1 \ with the remaining words providing cyclic permutations thereof. 
Again, the second trace is obtained by flipping signs in the exponents, and the difference of 
the two is 

- 2iT(l - T f sin V(t 3 - h) + 2iT 2 (l - T) sin V(t 3 - t x ) = -2iT(l - T)(l - 2T) sin V(t 3 - h) , 

(4.36) 
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plus cyclic permutations. Terms involving precisely one incoming currents yield II. Finally, 
terms containing more than one incoming current vanish. 

Each term of II exhibits two singularities not cancelled by the numerator. The terms 
differ however in regard as to whether the singularities may be attained within the region 
ti > ^2 > £3: both in the first term, neither in the second, and ti = t 2 + 21 but not t 2 = ^3 — 2/ 
in the last one. Within said region the regularization of the unaccessible singularities may be 
changed from — iO to +i0. The so modified expression II' is permutation symmetric, like I. 
This proves 

((T T11T12T13}} = 1 + 11' (4.37) 
for all ti, t 2 , t s > I. We can now compute its integral over < ti < t, {i = 1, 2, 3), for large t. 
The contribution from I is inferred from (I4.29p . The three terms of II' contain singularities 
in the variables t 3 , ti, respectively t 2 , which do not pinch the real axis. Hence II' does not 
contribute and the result is 

«(AQ,) 3 » = / dh «T JnW, 3 » = — T(l - T)(l - 2T) + o(t) . (4.38) 
Jo 27T 

4.4 Comparison of time orderings (continued) 

The purpose of this subsection is to show that eq. (13. 18f) persists upon replacing the currents 
by their second quantized counterparts: 

T* (Ife) • • • T(t k )) = f (/ft) ■ • • lift)) . (4.39) 

Inspection of the proof of Proposition [3] shows that we need to establish [Q(s), I+(t)] = 0, 
(t < s) and a similar property for i_ (t) . By time covariance it suffices to prove the claim for 
s = 0, where Q(0) = Q. Since by (I3.2ip we have 

[Q(s),J + (t)] = 0, (t<s), (4.40) 

at the level of first quantization, we have to make sure that the property is not destroyed by 
Schwinger terms: 

tr(pQ(s)p'I + (t)p) - tr(p'Q(s)pI + (t)p>) = , (4.41) 

which for s = indeed follows from pQp' = 0. Alternatively, one can verify that the Schwinger 
term does not change under a common translation of t and s, i.e. under conjugation of 
Q(s) and I+(t) by the propagator U(t). Schwinger terms satisfy the Hochschild condition 
S(AB, C) + S(BC, A) + S(CA, B) = and hence 

S([A, B},C) + S([B,C},A) + S([C,A},B) = 0. (4.42) 

As a result of (14.401) we have 

-i^S(Q(s+T)J + (t+T))=S([H,Q(s+T)]J + (t+T))+S(Q(s+T) (4.43) 
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